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Abstract. In this paper, the explicit form of maximal elements, known as 
shorted operators, in a subring of a von Neumann regular ring has been ob- 
tained. As an application of the main theorem, the unique shorted operator (of 
electrical circuits) which was introduced by Anderson- Trapp has been derived. 



1. Introduction 

Various partial orders on an abstract ring or on the ring of matrices over the 
real and complex numbers have been introduced by several authors either as an 
abstract study of questions in algebra, or for the study of problems in engineering 
and statistics (See, e.g. [T], [2J, 0], [7], [12], and [IS])- Also, a partial order on 
semigroups is studied by several authors (See, e.g. [5], [15], and [H]). In this 
paper we study the well-known minus partial order on a von Neumann regular ring 
which is simply a generalization of a partial order on the set of idempotents in 
a ring introduced by Kaplansky. For any two elements a, b in a von Neumann 
regular ring i?, we say a <~ b (and read it as a is less than or equal to b under the 
minus partial order) if there exists an x € R such that ax — bx and xa — xb where 
axa = a. Furthermore, we define the partial order <® by saying that a <® b 
if bR = aR © (b — a)R, and call it the direct sum partial order. The Loewner 
partial order on the set of positive semidefinite matrices S is defined by saying that 
for a, b 6 S, a <l b if 6 — a G S. The direct sum partial order is shown to be 
equivalent to the minus partial order on a von Neumann regular ring. It is known 
that the minus partial order on the subset of positive semidefinite matrices in the 
matrix ring over the field of complex numbers implies the Loewner partial order. 
The main result of this paper gives an explicit description of maximal elements in 
a subring under minus partial order (Theorem 13). As a special case, we obtain a 
result similar to the one obtained by Mitra-Puri ([IB], Theorem 2.1) for the unique 
shorted operator; which, in turn, is equivalent to the formula of Anderson- Trapp 
([2J, Theorem 1) for computing the shorted operator of a shorted electrical circuit 
(Theorem 17). 

2. Definitions 

Throughout this paper, R is a ring with identity. An element a S R is called von 
Neumann regular if axa = a for some x € R and x is called a von Neumann inverse 
of a. We will denote an arbitrary von Neumann inverse of a by a^\ An element 
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a G R is called weakly regular if xax = x for some x 6 R and x is called a weak 
von Neumann inverse of a. We will denote a weak von Neumann inverse of a by 
a,( 2 \ If ass = a and a;ax = x, then x is called a strong von Neumann inverse of a. 
We will denote a strong von Neumann inverse of a by a^ 1:2 \ A ring R is called von 
Neumann regular if every element in R is von Neumann regular. For convenience, 
we will use the terminology regular ring in place of von Neumann regular ring. For 
details on regular ring, the reader is referred to [5]. 

Let S be the set of all regular elements in any ring R. For a, b G S we say 
that a <~ b if there exists a von Neumann inverse a; of a such that ax — bx and 
xa = xb. This is known as the minus partial order as stated above for regular rings. 
The minus partial order clearly generalizes the definition of Kaplansky according 
to which if e, / are idempotents then e < / if ef = e = fe. 

We remark that for the ring of matrices over a field, it is known that a < 6 if 
and only if rank(b — a) = rank(b) — rank(a). 

Let T be a ring with involution *. If a; is a strong von Neumann inverse of a 
such that (ax)* = ax, (xa)* = xa and ax = xa then x is called the Moore-Penrose 
inverse of a and is denoted by a^ . Let M be the set of positive semidcfinitc 
matrices. For w G M and b G T, x is called the unique w-weighted Moore-Penrose 
inverse of b if x is a strong von Neumann inverse of b and satisfies (wbx)* — wbx and 
(wxb)* — wxb. For details on Moore-Penrose inverse, one may refer to Rao-Mitra 
[17] or Ben- Israel and Greville [3]- 

3. Preliminary Results 

The following result of Jain and Prasad ([8j, Theorem 1) will prove to be useful 
throughout this paper and, specifically, for providing an equivalent definition of the 
minus partial order on a regular ring. 

Theorem 1. Let R be a ring and let a, b G R such that a + b is a regular element. 
Then the following are equivalent: 

(1) aR(BbR = (a + b)R; 

(2) Ra © Rb = R(a + b); 

(3) aRnbR= (0) = Ra n Rb. 

From Rao-Mitra ([IT], Theorem 2.4.1, page 26), we have the following nice char- 
acterization of {c^ 1 -*} and {a^ 1,2 -*}. 

Lemma 2. Let R be a ring and let a G R. If x G {a*- 1 ^} ^ien {a^ 1 -*} = x + (1 — 
xa)i? + i?(l - ax). In addition, {a^ 1 * 2 )} = {a^W 1 )}. 

We now investigate properties of the direct sum partial order and its relation to 
the minus partial order. 

Let R be a regular ring. Recall a <® b if and only if bR = aR © (b — a)R. By 
Theorem 1, this is equivalent to Rb = Ra © R(b — a). It is straightforward to see 
that <® is a partial order. 

Next we show that the minus partial order is equivalent to the direct sum partial 
order on a regular ring. Hartwig-Luh showed that, when R is a regular ring, (2) is 
equivalent to (3) with the additional hypothesis that a G bRb (see [M], page 5). 

Lemma 3. Let R be a regular ring and a,b G R. Then the following are equivalent: 
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(1) a <® 6; 

(2) a <- 6; 

(3) {&«}C{a«}. 

Proo/. (1) =>■ (2) : As a <® b, bR = aR ® (b - a)R. It follows that aR C bR. 
Hence, a £ bR and thus a — bx for some x E R. As i? is a regular ring, for 
any g E {o^}, bgb — b. Thus oga = bg(bx) — (bgb)x = bx — a. Now aga = 
bga— (b~ a)ga — a— (b~ a)ga. Thus a — aga = (b — a)ga. But ai?n (b — a)R = (0) 
and a — aga = (6 — a)ga G ai? n (6 — a)R. Hence a — aga = and (6 — a)ga = 0. 
Therefore aga = a = oga and hence {o^ 1 -*} C {a' 1 )}. Indeed, this demonstrates 
that (1) (3). Now choose x = gag. Then aira = a(gag)a = aga = a and 
x G {a- 1 '}. Now 6a; = (bga)g — ag as bga = a. Furthermore, ax = agag = ag as 
aga = a. Thus ax = bx. Now bg(b — a) = bgb — bga = (6 — a) and (b — a)g(b — a) = 
bg(b — a) — ag(b — a) = (b — a) — ag(b — a). Hence ag(b—a) = (b— a) — (b— a)g(b— a) G 
aRn(b-a)R={0). Thus (b - a) = {b - a)g(b - a) and ag{b - a) = 0. It follows 
that agb = aga — a. Now xb = (gag)b = g(agb) = ga and xa = gaga = ga. 
Therefore xb = xa. Thus ax = bx and xa = xb for some x G {o^ 1 -*} and it follows 
that a <~ 6. 

(2) => (3) : This is well-known. We prove it here for completeness. As a < b, 
there exists some x G {a*- 1 -*} such that ax = bx and xa = xb. It follows that 
a = axa = bxa = axb and for any y G {b^}, aya = (axb)y(bxa) = ax(byb)xa = 
axbxa — (axb)xa — axa — a. Thus {o^} C {a*- 1 - 1 }. 

(3) => (1) : Given that {&«} C {a«}, ab^a = a for any 6« G {o«}. By 
Lemmai {&•»} = g + (1 - gb)R + R(l - bg) for g G For each a; G {M 1 )} 
there exists some r±, 7*2 G R such that a* = g + (1 — go)?*i + r%{Y — og). Multiplying 
on the left and right by a yields axa = a \g + (1 — go)ri + 7*2(1 — og)] a. Hence 
a = axa = a [g + (1 — go)?*i + 7*2(1 — og)] a = aga + a(l — gb)r\a + 07*2(1 — bg)a = 
a + a(l — gb)r\a + 07*2(1 — bg)a. Thus a(l — gb)r\a + ar-z(l — bg)a = 0. As 
a(l — gb)r\a + 07*2(1 — bg)a = holds for all 7*1 and 7*2, we can take, in particular, 
7*2 = which gives a(l — gb)r\a = for all 7*1 and hence a(l — gb)Ra = (0). 
Similarly, by taking 7*1 = 0, we conclude ai?(l — bg)a = (0). Now (a(l — gb)R) = 
(a(l-gb)R)(a(l-gb)R) = (o(l - gb)Ra)((l - gb)R) = (0) ((1 - gb)R) = (0). 
Similarly (i?(l — bg)a) = (0). Since R is a regular ring, it has no nonzero nilpotent 
left or right ideal. Thus, a(l - gb)R = (0) and R{1 - bg)a = (0). As 1 G R, 
a(l — gb) = and (1 — bg)a = 0. Therefore, bga = a = agb. Now for any 
ti, t<i G R, at\ = (bga)t\ — b(gat\) G bR and (6 — a)£2 = 0^2 — 0^2 = bt2 — {bga) £2 = 
b(t 2 -gat 2 ) ebR. Hence, ai?+ (6- a)i? C bR . Thus ai?+ (b~a)R = bR. Now we 
want to show that aRD (b — a)R = (0). For some u,v G R, suppose au — (6 — 0)7; G 
aR n (6 — a)i?. Then au = agau = ag(b — a)u = ag6t> — agav — av — av = as 
a = ago. Thus aR n (6 - a)i? = (0) and so bR = aR®(b- a)R. Hence, a <® 6 as 
required. □ 

We also note that proving directly (2) =£> (1) requires a brief argument. 



The Corollary that follows shows, in particular, that the minus partial order de- 
fined on the set of idempotents is the same as the partial order defined by Kaplansky 
on idempotents (See e.g. Lam [9], page 323). 
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Corollary 4. Let R be a regular ring and a,b 6 R such that b = b 2 . Then the 
following are equivalent: 

(1) o<- b- 

(2) a = a 2 = ab = ba. 

Proof. The proof is straightforward. □ 

Corollary 5. Let R be a regular ring and let a,b,c £ R with b = a + c. Then the 
following statements are equivalent: 

(1) a <~ b; 

(2) aRC\cR= (0) = Ra n Rc. 

Proof. It follows from Lemma 3 and observing that, in a regular ring, a <~ a + c 
if and only if a <® a + c if and only if (a + c)R = aR © cR. □ 



Hartwig ( 6 , Pages 12-13) posed the following questions, among others: 

(1) If R is a regular ring and aR fl cR — (0) = i?a fl Rc, does there exist a*- 1 -* 
such that a^c = = c a W? 

(2) Does a <~ c, 6 <~ c, aRCicR = (0) = RaHRc imply a + 6 <~ c? 

As a byproduct of the development of the direct sum partial order, we give 
an application that answers the above two questions of Hartwig. We do not know 
whether or not someone has answered these questions, as we could not find this 
in the literature. In any case, we believe that the answers we have given would 
be of interest to the reader. Below, we answer Question 1 in the affirmative and 
Question 2 in the negative by providing a counterexample. 

Proposition 6. (Hartwig Question 1) If R is a regular ring and aRC\cR = (0) = 
Ra fl Rc, for some nonzero elements a,c £ R, then there exists a nonzero a^ such 
that aWc = = caW . 

Proof. Let b = a + c. By Corollary 5, a <~ b. Then, by the definition of the 
minus partial order, for some a^\ ooP^ = ba^ and a^a = a^b. Now substituting 
b = a + c yields aa^ = (a + c)a (1) and a (1) a = a (1) (a + c). Thus aa (1) = aa^ +ca^ 
and a^a = a^'a + aS^c. It follows that ca^ = = a^c as required. □ 

Example 7. (Hartwig Question 2) 

Using matrix units eij, let a = en, b = ei^, and c — ei3 + en + e24. Clearly 
a <~ c and b <~ c. It is obvious that aR (~l bR — (0) = Ra n Rb. Since rank(c) — 
rank(a + b) = 2 — 2 = and rank(c — (a + b)) = 1, it follows that a + b c. 
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4. Main Results 

Let R be a regular ring and S be a subset of R. We define a maximal element 
in C = {x € S : x <® a} as an element b ^ a such that b <® a and if b <® c <® a 
then c = 6 or c = a. 

For fixed elements a,b,c G i?, we give a complete description of the maximal 
elements in the subring S = eRf, where e and / are idempotents given by eR = 
aR n cR and Rf = Ra (1 Kb. Here, C = {s 6 eRf : s <® a}. In the literature, 
maximal elements in C have been called shorted operators of a ([T], [5] and |13]V 

We begin with a result that is used frequently in the sequel. This is indeed con- 
tained in ( 15 1, Lemma 1) where the author proves the equivalence of 11 statements. 
However, for the sake of completeness, we provide a direct argument. 

Lemma 8. Suppose R is a regular ring and a, b e R such that {a«} n ^ 0. 

Then the following are equivalent: 

(1) aR C bR and Ra C Rb; 

(2) a <® b. 

Proof. Suppose aR C bR and Ra C -Rfo. It follows that a — rb — bs for some r, s £ 
i?. We claim that ab^a is invariant under any choice of b^\ Let x, y € 
be arbitrary. Now axa = (rb)x(bs) = r(bxb)s — rbs as 6x& = 6. Similarly, aya = 
(rb)y(bs) = r(byb)s — rbs as byb = b. Thus axa = aya for every x, y € {&^}- 
Hence ab^a is invariant under any choice of frW. Since we have assumed that 
{a (1) }n{6 (1) } ^ 0, there exists some g e {a (1) }n{o (1) }. Therefore a6 (1) a = aga = a 
for all &W. Hence C and by Lemma H a <® b. 

Conversely, if a <® b, then aR C oi? and Ra C i?6 follow by definition. □ 

We now demonstrate an important relationship between weak von Neumann 
inverses and strong von Neumann inverses under the direct sum partial order. 

Lemma 9. Let a £ R where R is a regular ring. Then the following are equivalent: 

(1) b is a weak von Neumann inverse of a; 

(2) There exists a strong von Neumann inverse c of a such that b <® c. 

Proof. Suppose 6 is a weak von Neumann inverse of a. For any fixed define 
u = at- 1 ' (a — aba)a^ and c = b + u. Then aca = aba + aua = aba + aa^aa^a — 
aa^abaa^a = aba+a — aba = a and cac = (b+u)a(b+u) = bab+bau+uab+uau = 
b + ba(aWaaW - a^abaa^) + (a (1) aaW - a^abaa^)ab+ 

(a^W 1 ) - a^abaa^)a(a^aa^ - a^abaa^) = b + baa^ - baa^ + a^ab - 
a^ab + a^aa^ -a^abaa^ -a^abaa^ +a^abaa^ = b + a^(a - a&a)a« = 
b + u = c. This shows that c is a strong von Neumann inverse of a. 

Now we want to show that b <® c. In other words, we will prove that bR(BuR = 
cR. Observe that cab = [b + a^(a — aba)a^]ab = bab + a^\ab — abab) = bab = b. 
Therefore b G cR. As c = b + u, it is clear that cR C bR + uR. As u — c — b 
and b € cR, uR C cR. It follows that cR = bR + uR. Now we want to show that 
bR n uR = (0). Let bp = uq £ bR n uR for some p,q € R. Multiplying ba on 
both sides yields bp = babp = bauq = ba[a^(a — aba)a^']q = (ba — baba)a^ q = 
(ba - ba)aWq = 0. Therefore bR PI uR = 0. Thus bR © uR — cR and we have 
demonstrated that b <® c. 

Conversely, suppose that there exists a strong von Neumann inverse c of a such 
that b <® c. As c is a weak von Neumann inverse of a, cac = c and thus a S {c*- 1 -*}. 



6 BRIAN BLACKWOOD, S. K. JAIN, K. M. PRASAD, AND ASHISH K. SRIVASTAVA 



By assumption b <® c and it follows from Lemma [3] that {c^} C {M 1 )}. Thus 
a G {c*- 1 '} C {b^ 1 '} and it follows that bab = b. Hence b is a weak von Neumann 
inverse of a. □ 

Lemma 10. Suppose R is a regular ring. Let y be a weak von Neumann inverse 
and z be a strong von Neumann inverse of an element a in the subring f Re such 
that y <® z. Then eyf <® ezf . 

Proof. Let a = fxe G fRe. Since y <® z, yR C zR and Ry C Rz. Thus, 
y — rz — zs for some r, s G R. It is straightforward to verify that zay = y = yaz. 
This gives (ezf)x(eyf) = (ezf)x(e(zs)f) = ez{fxe)zsf = ezsf = eyf. Similarly 
(eyf)x(ezf) — eyf. Thus (eyf)R C (ezf)R and R{eyf) C R(ezf). Asa = /xe 
is a common von Neumann inverse of y and z, it follows that (eyf)x(eyf) = eyf 
and (ezf)x(ezf) — ezf and so i is a common von Neumann inverse of eyf and 
ezf. By Lemma O ey/ <® ezf . □ 

Next, we give two key lemmas. We will assume throughout that a ^ S. 

Lemma 11. Let R be a regular ring. Then d G C is a maximal element in C if 
and only if for any d <® a such that dR C d R C eR, Rd C i?d C i?/, we /icwe 
d = df. 

Proof. Let d be a maximal element in C. If d is any element in R such that 
d' <® a and di? C d' R C e i?, i?d C i?d' C i?/, then clearly d' G ei?/. As d' <® a, 
d' G C. Then {a^} C {d^ 1 )} n {(d') (1) }. Hence, d <® d' by Lemmai Then by 
the maximality of d in C, d = d . 

The converse is obvious. □ 

Lemma 12. C = {eu/ : u is a weak von Neumann inverse of fa^e}. 

Proof. Let s = etf G C for some t G R. Then s <® a. By Lemma El {a^} C 
{ S W}. Therefore, we have (etf)a^ {etf) = (etf). In other words, (etf){fa^e)(etf) = 
(etf), proving that s = etf is a weak von Neumann inverse of fa^e. This shows 
that s = euf for some weak von Neumann inverse u of faP^e. 

Conversely, consider any u G (fa^e)^ and let x — euf. We want to show 
that x <® a. Now xa^x — (euf)a^ (euf) = eu (joS^e) uf — euf — x as 
ue(fa^e)^ 2 \ Hence {a (1) } C {a:^ 1 ^}. By LemmaEl x <® a and so x — euf G C. 

□ 

Theorem 13. maxC = {evf : v is a strong von Neumann inverse of fa^'e}. 

Proof. Suppose x — euf G C where u = (fa^ej . By Lemma [9l there is a 
strong von Neumann inverse v G eRf of foP^e such that u <® v and consequently, 
by Lemma 10, euf <® evf. Thus, we have x <® evf. Next, we will show that 
evf <® a. We have (evf)a^ (evf) = ev (fa^e) vf — evf as v G (/aWe)' 1,J l 
Hence {a^ 1 '} C {(eti/)^ 1 ^}. By Lemma [3 ew/ <® a. Thus maxC C {evf : v 
is a strong von Neumann inverse of/a«e}. Clearly, maxC is non-empty unless 
evf = a for each choice of v but this is not possible as we have assumed a $ S. 

Now suppose evf,ev'f G C such that v,v' are strong von Neumann inverses of 
fa^'e and evf <® ev' f. Therefore ev' fR — evfR © (ev' f — evf )R. Now we 
want to show that ev' fR = evfR. As evf, ev' f G C, evf <® a and ev' f <® a. 
Thus {a^} C {( e «/) (1) } and {a^} C {(ei;7) (1) }. So let a (1) be a common von 
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Neumann inverse of evf and ev' f. By assumption evfR C ev'fR. As shown in 
LemmaHni (ev' '/) a« (evf) = evf and (ev' f) (ev 1 f) = (ev'f). Now (ev'f) R = 
ev'fa {1) R = ev'fa^eR = ev' (f ev f a { ^ e) R C ev' f ev f R = evfR C eu'/K. 
Thus ev' fR = evfR. Similarly we can show that Rev' f — Revf. 

As Rev' f — Revf, we claim that ei// = evf. Let et// = revf for some r £ R. 
Now eu/ = ev'fa^evf — (revf)a^'evf = r(evf) — ev'f. Thus evf — ev'f. 
Hence maxC = {evf : v is a strong von Neumann inverse of/a«e}. □ 



We now provide an example to illustrate the previous theorem. 
Example 14. Note that we are choosing f to be of rank two. 
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Then one choice for a*- 1 -* 



For our choice of a 



strong von Neumann inverse of foS^e, we first choose its Moore-Penrose inverse 
and later its group inverse, as both are also strong von Neumann inverses. Let 
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Now ev\f < a because rank(a — evif) 
Thus evif SmaxC. 



We now find another element of maxC. The group-inverse v-i of faP~'e is 



V2 



Now ev2f < a because 



Iff II 00 

— - - --00 

8 8 8 o ■ | | o o 

lj |_0 1 

— ev2f) = 2 = 4 — 2 = rank (a) — rank(ev2f). Thus ev2f GmaxC. 



5. An Application 

In this section, as an application of our main theorem on maximal elements, we 
derive the unique shorted operator as of Anderson- Trapp (See [2], Theorem 1) that 
was also studied by Mitra-Puri (See [13], Theorem 2.1). We believe that there will 
be other such applications. 



8 BRIAN BLACKWOOD, S. K. JAIN, K. M. PRASAD, AND ASHISH K. SRIVASTAVA 



Throughout this section R will denote the ring of n x n matrices over the field 
of complex numbers, C. For any matrix or vector u, u* will denote the conju- 
gate transpose of u. In this section S will denote the set of positive semidefinite 
matrices. 

Recall, the Loewner order, <l, on the set S of positive semidefinite matrices in 
R is defined as follows: for a, b £ S, a <l b if b — a £ S. 

Suppose a £ S and c £ R. As in the previous section, eR = aR fl cR, e = e 2 , 
and choose f — e* . Clearly, / £ Ra because a is hermitian. Let Cl — {s £ 
eRf fl S : s <l a} = {s £ e5/ : s <£ a}. 

Under this terminology, the set C in the previous section will become, C = {s £ 
eS7 : s <® a}. 

We will assume that a ^ eSf. This is equivalent to the assumption that 
rank(e) ^ rank(a), as shown in the remark below. 

Remark 15. rankie) — rank(a) if and only if a £ eSf. 

Proof. Suppose rank{e) = rank(a). So eR = aR as eR C aR. Then a = ex for 
some x E R and by taking conjugates, a = x*e* , i.e., a € i?e*. Hence, a £ eRe*. 
As a g 5, a £ S fl ei?e* = eSe*. For if exe* £ 5 then exe* = e (exe*) e* £ e5e* 
and so S fl ei?e* C eSe*. The reverse inclusion is obvious. 

Conversely, suppose a £ eS 1 /. As eR = aR fl cR, we have e = ax and so 
rank{e) < rank(a). As a 6 eS*/, a = ese* for some s £ S. Therefore rank(a) < 
rank(e). Hence, rank(e) = rank(a). □ 

The following lemma is folklore. 

Lemma 16. Suppose a,b £ S. If a <® b then a <£ fo. 

Proof. Suppose a <® b. Equivalently, (b — a) <® b and by Lemma[3]we know that 
{&«} C {(6-a) (1) }. Thus, fet 

is a von Neumann inverse of (b — a). From 
[llj . as 6 is positive semidefinite, 6^ is positive semidefinite. Thus b — a = 
{b - a) {b - a) > L 0. Hence (b - a) E S and a < L b. □ 

Theorem 17. Let a E S and let /J &e i/ie a-weighted Moore-Penrose inverse of 
f. Then maxC = maxCi = {a/d/}- 

Proof. By Theorem 13, maxC = {evf : v is a strong von Neumann inverse of 
/ct^e}. By assumption, e E ai? and so e = ax for some x £ R. By taking 
conjugates, e* = x*a as a <E S. In addition, as / £ i?a, f — ya for some y <E R. 
This yields that fa^e — yaa^ax — yax and thus fa^e is independent of the 
choice of a^ 1 - 1 . We may then choose the Moore-Penrose inverse for a^ . Next, we 
want to show that a strong von Neumann inverse of fa^e is also unique. Note that 
fa^e — e*a i e is positive semidefinite, as the Moore-Penrose inverse of a positive 
semidefinite element is positive semidefinite [TT]. As a £ S, we can write a = zz* 
for some z £ R. Now fR = yaR = yaa^aR = fa^aR = fa^R = fzz*R = 
fzR = (fz) (fz)* R = fzz*f*R = fa^eR. Similarly Re = Rfrfe. It follows 
that / = fa'ep and e = qfa^e for some p,q £ R. Consider an element evf £ 
maxC. Then evf = qfa^evfa^ep = qfa^ep, showing that evf is independent of 
the choice of strong von Neumann inverse v of faJe. Thus maxC is a singleton set 

consisting of the element e (fa^e) f . Since a £ S, a) £ S and hence e (e*a i 'e) / = 

e(/at e ) t /£^. 
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Next, we proceed to show that maxC = {aflf} also. Recall that aflf is 
hcrmitian and so aflf = (aflf)* = f* (fl)* a* = (/]/)* a. Since f\f is an 
idempotent, we get aflf = o(/t/)(/t/) = (/+/)* o(/+/) and thus aflf G 5. 

We now prove that a/2/ <® a. Let a*- 1 -* be an arbitrary von Neumann inverse 
of a. Then (aflf) a^ (aflf) - (afl)(ya) a W (aflf) = {afly)aa^a(flf) = 
aflyapj = aftfftf = aflf. Hence {a«} C {(a/t/) (1) }. Consequently, by 
Lemma 3, aflf <® a which gives aflf G C. 

Furthermore, by Lemma 16, aflf <® a gives aflf <l a and hence aflf G Cl- 

Finally, we show that for every d G Cl, d <l aflf- As d G S C i?/, write d = u/ 
for some u E R. Then d/j/ = ufflf = uf = d = (flf)*d (flf) as d is hermitian. 
Now consider a/]/ - d - (/]/)* a (/]/) - (/]/)* d (/]/) = (/]/)* ( a - d) (/]/) , 
which is positive semidefmite and thus aflf — d G S. Hence d <l aflf. 

Thus aflf is the unique maximal element in Cl provided aflf ^ a. We have 
shown above that aflf G Cl and thus aflf G eSf. But by assumption a £ eSf . 
So aflf ^ a. Therefore, aflf is unique maximal element in Cl and it also belongs 
to C as we have already proven that aflf <® a. 

Now, because e (/c^e) / is the unique maximal element in C and aflf G C, 
o/J/ <® e (/at e ) f / . By Lemma 16, aflf < L e (/ a t e ) t / as e (/a^ / G Cl- 
We have shown above that for every element d G Cl, d <l aflf and thus a/J/ = 
e(/at e ) f /. Hence, maxC = max Cl = {»/!/} as desired. □ 



The following examples demonstrate the result proved in the previous theorem, 
i.e. aflf = e (fa^e) f and so maxC = max Cl = {aflf}. Furthermore, maxC 
agrees with the formula given by Anderson- Trapp for computing the shorted oper- 
ator as when we are given the impedance matrix a. 

The Anderson- Trapp formula states that if a is the nxn impedance matrix then 
the shorted operator of a with respect to the fc-dimensional subspace S (shorting 

an — ai2c4 2 «2i 




k ports) is given by as 



where a is partitioned as 



such that an is a k x k matrix. We show that the maximum 



an ai2 

a = 

^21 Cl22_ 

element aflf obtained by us is permutation equivalent to as, i.e. P T aflf P = as 
for some permutation matrix P. 



Example 18. Let e 



i 

2 



1 

2 













and then f = e* 















1 



1 
1 



Then one may check that fl = f 



i 

2 


1 

2 


1 

2 


1 

2 





So aflf 



10 10 



10 10 

1 











1 

2 u 



h o 

1 



Suppose 



We now show that aflf = e (fa^e)^ f. Now, a* 
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1 

4 





1 

4 








1 








1 

4 





1 

4 














1 



and (/a^e)^ = 



I— 1 





1 

















I— 1 





1 














1 



Thus e (f 



1 





1 

















1 





1 














1 



Hence e (/a^e) 1 ^ / 



1 





I— 1 

















1 





1 














1 



a flf as proved in the theorem. We may 



verify that aflf <® a. This follows from rank(a) — rank(aflf) = 3 — 2=1 = 
rank(a — aflf). We know then aflf <l a. Thus maxC = maxCx = {aflf}. 
We now compute the shorted operator as given by Anderson-Trapp. We partition 



a as follows: a 



Then as 







1 





1 














1 














1 





1 











[o 





oj 




[1 


"1 





1" 




"0" 







1 










[1 


1 — 1 





1 


















"1 





1 













1 












1 





1 























Now for 



1 — 1 




















1 








1 








1 









PaflfP T = a s . 
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